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I. INTRODUCTION 



A challenging task in high precision calculations of energy levels of few electron atoms 
or ions is the accurate solution ip of the nonrelativistic Schrodinger equation. This wave 
function ip is next used to obtain relativistic and quantum electrodynamics (QED) effects 



including finite nuclear mass corrections, 
electron wave function ip achieved so far 1 

ip = ^Ci4>i 



he most accurate representation of the three- 
uses the Hylleraas basis set 3], namely 



-wi n-w 2 r 2 —w 3 r 3 m n 2 r n 3 n 4 n 5 n 6 
'23 '31 '12 ' 1 '2 '3 ) 



(1) 



with rii being nonnegative integers. This basis set allows the Kato cusp condition to be 
satisfied to a high degree, thus ensuring good convergence for matrix elements of relativistic 
and QED operators. All these matrix elements can be expressed in terms of the Hylleraas 
integrals /, see Eq. (j2j). The calculation of the nonrelativistic energy involves integrals with 
nonnegative n^. Such integrals have been worked out by King et al. in the series of works 
4], and more recently by Yan and Drake 1,5], Sims and Hegstrom 6], and by us in Ref. 
Our approach relies on analytic recursion relations, which are highly efficient and sufficiently 
stable to achieve accurate numerical results. 

Matrix elements of relativistic operators involve the so called extended Hylleraas integral, 
where one of n 7 - = —1 (see Eq. (12"D). These integrals have been worked out first by King et 

n u nn 

al. in |8j, and later by Yan and Drake in Refs. [9( and by us in [2j, [10[]. The use of recursion 
relations allows for a significant increase in the size of the basis set and the accuracy of 
the obtained results |2|. Finally, matrix elements of QED operators involve integrals with 
rti = —2. Such integrals have been worked out by Yan in but no numerical results 
for any particular integral have been published, which would serve as a test of achieved 
accuracy. In this work we develop recursion relations for three-electron Hylleraas integrals 
either with \jr\- or with 1/rf, and obtain a one dimensional integral representation for the 
master integral which is suitable for precise numerical evaluation. 

In Section II we recall known results for the regular Hylleraas integrals, in Section III we 
treat integrals involving and in Section IV we treat integrals involving 1/r?-. Apart 

from known results, we present also a new form of the master integral with the hope that it 
can be used in the future to obtain a fully analytic result. In Section V we develop recursion 



2 



relations for integrals involving 1/rf which are very similar to those involving 1/rf. In 
Section VI, which is the most difficult one, we obtain a full set of recursions for Hylleraas 
integrals with 1/r?-, and in Section VII we present numerical results together with a short 
summary. 



II. REGULAR THREE-ELECTRON HYLLERAAS INTEGRAL 



The regular three-electron Hylleraas integral is 
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with all rij > 0. The most convenient way to perform this integral is by using recursion 
relations. The initial values with rii — 0,1 are known explicitly [7|, Il2l |. 
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where 



Lg(z) = Li 2 (l -x) + Li 2 (-x) + ln(x) ln(l + x) , (11) 

and Li 2 (x) is the dilogarithmic function. All other Hylleraas integrals for arbitrary but 
nonnegative n, can be obtained by 6 independent recursion relations for each rij, see Ref. 

The numerical stability of these recursions is an issue. For the maximal number Q = 
J2i n<i = 30 we used sextuple precision arithmetic and verified against octuple precision that 
results for nonrelativistic energies are numerically significant for at least 16 digits. 



III. EXTENDED HYLERAAS INTEGRAL WITH l/rf 



More difficult to evaluate are extended Hylleraas three-electron integrals involving addi- 
tional single powers of or r^- in the denominator. The first kind of integral with l/rf can 
be obtained by integration with respect to a corresponding parameter u>j, namely 

POO 

f(0, 0, 0; -1, n 5 , Tie) = / dw, f(0, 0, 0; 0, n 5 , n 6 ), (12) 

where /(0, 0, 0; 0, 715, r^) is obtained by recursions which are numerically stable at hig h values 
of w\. The integration in Eq. ffl2l) is per formed using the adapted quadrature [13| for the 
class of functions involving logarithms }13| . which allows one to achieve high precision at low 
evaluation cost. The function /(rij, n 2 , n 3 ; — l,n 5 ,n 6 ) for arbitrary integer values of ni,n 2 
and n 3 can be obtained Q by recursion relations starting from /(0, 0, 0; —1, n 5 , n 6 ). 



IV. EXTENDED HYLERAAS INTEGRAL WITH 1/r 2 



The second kind of extended Hylleraas integral involves 1/r?-, and we will investigate it 
here in more detail. Recursion relations start from the master integral 



h(w u w 2 ,w 3 ) = /(-l, 0,0; 0,0,0) 

d 3 T\ f d 3 r 2 f d 3 r 3 
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and some other simpler two electron-like integrals. By using integration by parts identities 
[ijj], the following differential equations for the h function have been obtained in Ref. 10] 
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It is convenient to express the function h 



h(w 1 ,w 2 ,w 3 ) 



1 
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in terms of a dimensionless function H(x,y) of two variables x and y. Then the first two 
differential equations become apparently equivalent and take the form 



dH 1 / 1\-V2 



where 



F = F(w 1 ,w 2 ,w 3 ) = F[j4 : x + y + ^,y/y,-^= 



and the third differential equation, with respect to Wi becomes 
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where 



R'(x,y) = ^ W2 3 Wi R{w 1 ,w 2 ,w 3 ) 



Ax + y + ^j ' R^4x + y + ^,y/y,^= 



(24) 



The homogenous differential equation in fl23|) is satisfied by complete elliptic integrals 
K(l/2 ±x), and the solution of the inhomogeneous equation is obtained by Euler's method 
of variation of constants [lo| 



H(x,y) = - 

71 



+K {\- x )L /2 dzR, ^ K {\ +z 
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This integral form is convenient for the numerical calculation of the master integral 
/i(wi,w 2 ,w 3 ) and its derivatives. However, a different and apparently simpler form can 
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be obtained from Eq. (T2TT) 

= r — F^4x>+(u-l/uy,l/u,u), (26) 

Jo u ^4x' + (u- 1/uf 

where x' = x + 1/2. This integral can be further transformed to a form involving Jacobi 
elliptic functions, but we have not been able to perform it analytically. In particular cases, 
where x = ±1/2, this integral becomes (w 2 < W3) 



■1 ^ = r^ d i F{1+<iy)i (27) 



2 ' W3 J Jo 1 + u 



H(~,-)= r du^— 2 F(l-u\l,u\ (28) 
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-, — ) = du , 

2 W3J Jo 1 — 

and this can be expressed in terms of Li 3 , Li 2 , and logarithmic functions. We use Eqs. (I27|l28p 
in cases w\ > W2 + and w\ < \w2 — w^\, and for its numerical calculations we employ 



adapted quadrature 



13J with 120 points, to achieve about 64 significant digits accuracy. 



Recursion relations for /(— 1, n 2 , 713, n^, n 5 , tiq) have been derived in [lOj. Among them, 
the most numerically unstable is the one which increases the parameter 77,4, see Eq. (36) of 
Ref. [lo|. Namely, close to singularity points, x = ±1/2 and for small u>i, we use Taylor 
expansion to avoid numerical instabilities. 



V. SINGULAR HYLLERAAS INTEGRAL WITH 1/rf 



The singular Hylleraas integrals which involve 1/rf and 1/rf- are needed for the compu 



tat ion of QED effects 



111 ]. Let us first define a distribution P(l/r 3 



lim J d 3 r (j)*{r) 
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Any factor 1/r 3 in the following will be understood to be defined in the above sense, and 
we will drop the symbol P(). It follows from this definition that 

= C \n(w) + 0{l/ u) (30) 



for large uj. Moreover, we will use in the derivation below the following integral representa- 
tions 
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To obtain recursion relations for the singular three-electron Hylleraas integral one first 
considers the integral G 

G(mi,m 2 ,m 3 ;m4,m5,rne) = —— f d 3 k x [ d 3 k 2 f d 3 k 3 (k\ + u\y mi (k\ + u\y m2 
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which is related to / by /(0, 0, 0, 0, 0, 0) = G(l, 1, 1, 1, 1, l)| ul=U2=1i3= o- The following 9 
integration by parts (IBP) identities are valid because the integral of the derivative of a 
function vanishing at infinity vanishes, 
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= id(i, j) = J d 3 h J d 3 k 2 J d 3 k 3 [kj (hi + u{) 
(kj + ujy 1 (kl + ul)-\kl 2 + w{Y l (kl 3 + wlY 1 (k\ x + wl)-'] , (33) 



where i, j = 1, 2, 3. The reduction of the scalar products from the numerator leads to identi- 
ties for linear combination of the G functions. If any of the arguments is equal to 0, then G 
becomes a known two-electron Hylleraas type integral. These identities can be used to derive 
various recursion relations. Here we derive a set of recursions for /(ni, n 2 , TI3, — 2, 715, no). 
This is achieved in a few steps. In the first step we use integration by parts identities in 
momentum representation Eq. fl33|) . to form the following linear combination 

id(2, 2) + id(3, 3) - id(l, 1) = 2 [G(0, 1, 1, 1, 1, 2) + G(0, 1,1,1, 2, 1) - G(l, 0, 1, 1, 1, 2) 

-G(l, 1, 0, 1, 2, 1) - G(l, 1, 1, 1, 1, l)/2 - G(2, 1, 1, 1, 1, 1) u\ 
-G(l, 1, 1, 1, 1, 2) (ul - u\) + G(l, 2, 1, 1, 1, 1) u\ 
-G(l, 1, 1, 1, 2, 1) (ul - ul) + G(l, 1, 2, 1, 1, 1) ul 
+G(l,l,l,2,l,lK] = 0. (34) 

We use Eq. ( 13~TI) to integrate with respect to w\ and differentiate over ui, w 2 , u%, w 2 , and 
W3 to obtain the main formula 
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where F is the two electron Hylleraas integral, defined in the appendix. The general formula 
in Eq. fl35|) does not work in the case 1 + n\ = n 2 + n 3 . In this special case we use IBP 
identities in coordinate space and limit ourselves only to identities of the form 
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The identities id(2) and id(3) 
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replace the main recursion in Eq. ( 1351) for the case 1 + ni = n 2 + n 3 , and can be used also 
for all other 7ij under the conditions that n\ > 0, n 3 > or ni > 0, n 2 > 0, respectively. 
The case of n\ = 0, n 2 = 0, n 3 = 1 or ni = 0, n 2 = 1, n 3 = 0, which is not covered by above 
recursions, remains to be treated. Therefore, in the third step one obtains the necessary 



recursions from the following combination of IBP identities in the momentum space: id(l, 1) 
+ id(2,l) + id(3,l) =0 

/(0, 1,0, -2, n 5 ,n 6 ) = — \ 
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This completes the evaluation of singular Hylleraas integrals involving 1/rf. 
VI. SINGULAR HYLLERAAS INTEGRAL WITH l/r|- 



The derivation of recursion relations and the master integral is similar to the one for 



Hylleraas integral with 1/rE, see Ref. [10|]. In the first step of deriving recursion relations 



we take the difference id(3, 2) — id(2, 2) and use it as an equation for G(l, 2, 1, 1, 1, 1) 
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Similarly, the difference id(2, 3) — id(3, 3) is used to obtain (7(1,1,2,1,1,1). These two 
equations are used now to derive recursions in n 2 and n 3 . With the help of Eq. (I3"T|) one 
integrates Eq. (j4"4"j) with respect to u\, which lowers the first argument ni to —2. Next, one 
differentiates with respect to u 2 , w 3 , Wi, w 2 , W3 at u 2 = u 3 = to generate arbitrary powers 
of r 13 , r 12 , r 1; r 2 , r 3 and obtains the recursion relation in n 2 
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-2 n 2 n 6 wi /(-2, n 2 - 1, n 3 , n 4 + 1, n 5 , n 6 - 1) 

+2n 2 wi w 3 f(-2,n 2 - l,n 3 ,n 4 + l,n 5 ,n e ) 
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+2n A wiw 3 f(-2,n 2 + l,n 3 ,n 4 - l,n 5 ,n 6 ) 

+n 6 wl f(-2, n 2 + 1, n 3 , n 4 , n 5 , n 6 - 1) 

~{n 2 + n 6 ) <5(n 4 ) T(n 3 + n 5 - 1, n 2 + n 6 - 1, -2, w 2 , u> 3 , 0) 

+u? 3 5(n 4 ) T(n 3 + n 5 - l,n 2 + n 6 , -2, u> 2 , w 3 ,0) 

+n 6 5(n 3 ) T(n 4 + n 5 - l,n 6 - l,n 2 - 2, u>i + w 2 ,w 3 , 0) 

-u> 3 <5(n 3 ) T(n 4 + n 5 - 1, n 6 , n 2 -2,wi + w 2 , w 3 , 0) 

+n 2 5(n 5 )T(n 6 -2,n 3 + n 4 -l,n 2 -l,w 3 ,w 1 ,0) , (45) 

where n« > 0, and the formula for /(— 2,n 2 ,l + n 3 ,n 4 ,n5,ne) can be obtained from the 
above one using symmetries of the function /. These recursions assumes that values of 
/(— 2, 0, 0; n 4 , n 5 , n e ) are known. We again obtain them using IBP identities. These are 9 
equations, which we solve against the following X i= i t9 unknowns at u 2 = u 3 = 

X 1 = G(l,2,l,l,l,lX 

X 2 = G(l, 1,2, 1,MK, 

X 3 = G(l, 1,1, 1,2,1) uj, 

X A = G(l,2,l,l,l,2)^, 

X 5 = G(l,2,l,l,l,l), 
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X 6 = G(l, 1,2, 1,1,1), 

X 7 = G(l,l,l,2,l,l), 

X 8 = G(l, 1,1, 1,2,1), 

X 9 = G(l,l,l,l,l,2). (46) 



Equations for X 7 , X$, and Xg are 



(3wl + w\ - wl) G(l, 1, 1, 1, 1, 1) -2u\w\ G(2, 1, 1, 1, 1, 1) 

-Aw\ ™l G(l, 1, 1, 1, 2, 1) - 2 w\ {w\ + wl - wl) G{1, 1, 1, 2, 1, 1) 

+F 1 (u 1 )w 2 - F 2 ( Ul )w 3 , (47) 

(3 wl - wl + wl) G{1, 1, 1, 1, 1, 1) - 2 u\ wl G(2, 1, 1,1,1, 1) 

-Awl wl G(l, 1, 1, 1, 1, 2) - 2 wl (wl - wl + wl) G(l, 1, 1,2,1, 1) 

-F 1 (u 1 )w 2 + F 2 {u 1 )w 3 , (48) 

(3 w\ - 4 wl wl + w\ - 4 wl w\ - 2 w 2 2 w 2 3 + wf) G(l, 1,1,1, 1, 1) 

-2 wl {w\ -2wlwl + wl-2wlwl-2w 2 2 wl + w\) G(l, 1, 1, 2, 1, 1) 

-2wl(ul (wl- wl- wl) + 2wlwl)G(2, 1,1, 1,1,1) 

-F^Ui) w 2 {wl -w 2 2 + wl) - F 2 (ui) w 3 {wl + wl- wl) , (49) 



where 



F 1 {u 1 ) = 2 w 2 [G{1, 1, 1, 0, 2, 1) - G{1, 1,1,1, 2, 0) + G{2, 0, 1, 1, 1, 1) - G{2, 1, 1,1,1, 0)] , 

(50) 

F 2 { Ul ) = 2 w 3 [G{1, 1, 1, 0, 1, 2) - Gil, 1,1,1, 0, 2) + G(2, 1, 0, 1, 1, 1) - G{2, 1, 1, 1, 0, 1)] . 

(51) 

One performs the u\ integration and obtains 

= F 1 w 2 - F 2 w 3 + {wl + wl - wl) f{-2, 0, 0, 0, 0, 0) 

-2 wl w 2 f{-2, 0, 0, 0, 1, 0) - Wl {wl + wl - wl) f{-2, 0, 0, 1, 0, 0) 
+wlT{0,-l,-l,w 1 ,w 2 + w 3 ,0), (52) 

= -F 1 w 2 + F 2 w 3 + {wl - wl + wl) f{-2, 0, 0, 0, 0, 0) 

-2 wl w 3 f{-2, 0, 0, 0, 0, 1) - Wl {wl - w 2 2 + wl) f{-2, 0, 0, 1, 0, 0) 
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+w\ T(0, -1, -1, w 1 ,w 2 + w 3 , 0) 



(53) 



-F 2 w 3 {w\ + w\ - w\) - F l w 2 {w\ - w\ + w\) + w f(-2, 0, 0, 0, 0, 0) 
-w Wi /(-2, 0, 0, 1, 0, 0) - 2 w\ w\ w\ f(0, 0, 0, 0, 0, 0) 
+wf (wf -w 2 2 - w l) r (°; -!> -!> ™i> w 2 + w 3 , 0) , 



(54) 



where w is defined in Eq. ([Til and 



/■oo 

F = / d«i« a -Fj(ui) 
Jo 



in the sense of the integral defined in Eq. ( 13TI) , with the result 



Fi 



lnfl + ^ lnf^-ln 2 ^)+2 1n 



(55) 



In 



-Li 2 



U>2 



Wi + w 3 
w 2 



In 1 



w 2 



Wi + w 3 



2Li 2 ( 

Wi 



-2 ln(w 3 ) lnf Wl 
W 2 +W 3 J \w 2 + w 3 

w 2 



Li 2 |l-^)-Li 2 (-^ 



-Li 2 1 



Wi+W 3 J \ w 1 + w 3 , 

r(o, -l, -2, w 2 , ™ 3 , o) - r(o, -l, -2, o, w 2 ) 



-w 2 T(0, 



•1, 0, wi, w 2 ) + w 2 T(0, -1, -1, w 2 , w x + w 3 , 0) 



(56) 
(57) 



lnfl + ^ ln(^-ln 2 ^)+2 1n 

V Wo J \w 2 J 



In 



-Li 2 



w 3 



Wx + w 2 
w 3 



.w 2 , 
\yy( 1 



Wl 



Wi 



w 3 



Wx + w 2 



-2Li 2 



w 2 + w 3 

W3 
W X 



— 2 ln(w 2 ) In 

Li 2 ( 1 - — ) - Li 2 

w 2 



Wi 



w 2 + w 3 

W3 
W 2 



Li 2 1- 



w 3 



Wi+W 2 / \ Wi + W 2 , 

r(o, -l, -2, W3 , w 2 , o) - r(o, -l, -2, o, Wl , w 3 ) 

-w 3 r(o, -l, -l, o, Wi, w 3 ) + w 3 r(o, -l, -l, w 3 , Wi + w 2 , o) . 



(58) 
(59) 



Under the replacement /(— 2, 0, 0, 0, 0, 0) = h/wi = h(0, 0, 0)/wi one obtains simplified 

differential equations 

dh dh 
= 2wiW 2 h {wl + w\- w\) h Fi w 2 - F 2 w 3 + wl T(0, -1, -1, w u w 2 + w 3 , 0) , 

OW 2 ~ OWi 



= 2 Wl w 3 ^ + {wl-wl + w 2 ^ <)h '' " 2 



dh 
dw 3 



dwi 



(60) 

Fi w 2 + F 2 w 3 + wl T(0, -1, -1, Wi, w 2 + w 3 , 0) , 



dh 



= w - F 2 w 3 (wl + w 2 2 - wl) - Fi w 2 (wl -wj + wl)-2 wl w\ wl f(0, 0, 0, 0, 0, 0) 
+wl (wl -w\- wl) T(0, -1, -1, Wi, w 2 + w 3 , 0) . 



(61) 
') 

(62) 
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One differentiates Eq. (!60|) with respect to w±, w 2 , and w 3 , and obtains in this way recursions 
for h in 77,5 

h(n 4 ,n 5 + l,n 6 ) = -— ^ [-n 5 Fi(n 4 , n 5 - 1, n 6 ) + w 2 Fi(n 4 , n 5 , n 6 ) 

ZW1W2 L 

+n 6 F 2 (n 4 , n 5 , n 6 - 1) - u> 3 F 2 (n 4 , n 5 , n 6 ) 

+ (n 4 - 1) n 4 r(— 1, n 4 - 2, n 5 + n 6 - 1, 0, w 1: w 2 + w 3 ) 

-2n 4 w 1 T(-l,n 4 — l,n 5 + n 6 — l,0,w;i,u7 2 + w 3 ) 

+w\ r(— 1, n 4 , n 5 + n 6 — 1, 0, wx, w 2 + w 3 ) 

— n 4 (n 4 + 2 715 — 1) fo(n 4 — 1, n^, tiq) 

+2 n 4 u? 2 h(ri4 - 1, n 5 + 1, n 6 ) 

+2 (n 4 + n 5 ) wi /i(ri 4 , n 5 , n 6 ) 

-(n 5 - 1) n 5 h(n 4 + 1, n 5 - 2, n 6 ) 

+2 n 5 w 2 h(n 4 + 1, n 5 — 1, n 6 ) 

+ (n 6 - 1) n 6 /i(n 4 + 1, n.5, n 6 - 2) 

-2 n 6 u? 3 /i(n 4 + 1, n 5 , n 6 - 1) 

— (w\ + w\ — wl) h(ri4 + 1, n 5 ,n 6 ) , (63) 



where 



Fi(ni,n 5 ,n 6 ) = 5(n 4 ) T(-2, n 5 , n 6 - 1, 0, w 2 , w 3 ) 

-n 5 T(-l, n 5 - 1, n 4 + n 6 - 1, 0, u> 2 , wi + w 3 ) 

+w 2 n 5 , n 4 + n 6 - 1, 0, u> 2 , u>i + w 3 ) 

+(n 5 - 1) 5(n 6 ) T(0, n 4 - 1, n 5 - 2, 0, wx, w 2 ) , 

-w 2 5(n 6 ) r(0, n 4 - 1, n 5 - 1, 0, w h w 2 ), (64) 

F 2 (n 4 ,n 5 ,n 6 ) = 5(n 4 ) T(-2, n 5 - 1, n 6 , 0, w 2 , w 3 ) 

-n 6 T(-l, n 6 - 1, n 4 + n 5 - 1, 0, w 3 , u> 4 + w 2 ) 

+w 3 T(-l, n 6 , n 4 + n 5 - 1, 0, w 3 , w x + tu 2 ) 

+ (n 6 - 1) 5(n 5 ) T(0, n 4 - 1, n 6 - 2, 0, Wi, w 3 ) 

-w 3 5(n 5 ) r(0,ra 4 -l,n 6 - 1, 0, tox, ty 3 ) , (65) 
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and the formula for h(n&, TI5, 1 + ne) is analogous. Eq. (|62|) is transformed now to the more 
explicit form 

dh gi g 2 g 3 # 4 



+ 



+ 



+ 



dwi w\ + w 2 + w 3 w\ — W2 — w 3 wi — W2 + w 3 wi + W2 — w 3 



(66) 



9i 



92 



1 

4 

+Li 2 



7T 



+ In 

6 v lt'2 

w 2 



+ Li 2 



w 3 
w 3 



+ [ln(w 2 w 3 ) - 2] In 



W2 + W 3 



-lnf— ) ln(^ + [2-ln(w 2 w 3 )] lnf- 
Vw 2 / Vw 3 / V'' 



Wi 



Li; 



,w 2 
w 2 

Wl 



Li 2 



W 3 / ' "' V w 2 + w 3 , 

wi/ Vwi + w 2 / VW1+W3 



(67) 



(68) 



93 



-_ + 2Lg(— ^— 
6 Vwi + w 2 



-Li 2 - 



w 2 
w 3 



2Li 2 



W3 
Wi 



+ 2Lg 
hLi 2 



w 2 + w 3 

U>3 
W 2 



2Li 2 (-^)+Li 2 (;i-^ 

Wi/ V w 3 



Li 2 1- 



W3 
W 2 



(69) 



0.i 



6 Vwi + w 3 



+ 2Lg 



it>i 



w 2 + w 3 



2Li 2 (-^-Li 2 fl-^ 

wi/ V w 3 



+Li 2 - 



w 2 
w 3 



2Li 2 



w 3 



_ Li2 _^ + Li 2 (i-H 

Wi/ V w 2 / V w 2 



(70) 



These formulae are used to obtain forward recursions for h{n^, 0, 0) similar to those used 
in two-electron Hylleraas integrals Namely, if we denote by hi = gi/(w\ + . . .) the 

corresponding term in Eq. (166]) . then the recursion for example for h 2 (n) defined by 

d 



h 2 (n) 

in the general case is the following 

h 2 (n) 



h 



2 • 



dwi. 

g 2 (n) +nh 2 (n - 1) 



(71) 



(72) 



Wi - w 2 - w 3 

In the special case of Wi ~ w 2 + w 3 one calculates at first recursions exactly at wi = w 2 + w 3 

92(n + 1) 



h 2 (n) 



n + 1 



(73) 



and obtains h 2 (n) for wi 7^ w 2 + w 3 using the Taylor formula. For this one notes that first 
derivatives of <?, are quite simple, 
dgi 1 



<9wi 4 wi 



[ln(wi + w 2 ) + ln(wi + w 3 ) - 2], 



(74) 
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dg2_ 



dg 3 



[ln(wi + w 2 ) + \n(wi + w 3 ) - 2] + 



1 



In 



( — — — 

\Wi + Wo/ 



In 



W2 



4(wi-W 2 + W 3 ) 4 (u^ + U> 2 + U?3 



4(^i + U7 2 -W 3 ) \W 1 + W 2 

' W W3 Uln 



(75) 



Vwi + w 2 



w 2 



W\ + w 3 



In 



Wl + w 3 



In 



u> 3 



Awi \Wi + w 2 J A(wi + w 2 — w 3 ) \wi + w 2 



(76) 



1 



In 



4 — w 2 — W3) \w 2 + u;3/ 4 (wi + w 2 + U7 3 ) 



1 



In 



dg 4 



In 



»i + w 2 



In 



w 2 



dwi Aw\ xwi + w?,/ 4(wi — w 2 + w 3 ) \wi + w 3 

Wi 



In 



A(w 1 — W2 — w 3 ) \w 2 + w 3 J A(wi + w 2 + w 3 ) 



+ 



1 



In 



^3 
W\ + w 2 



w 2 
W\ + w 3 



+ In 



Wi 



w 2 + w 3 



(77) 



+ In 



w 2 + w 3 



Independly of recursions one needs the function h. It is obtained from Eq. ([66 

dh 



dw\ 



dw\ 



(78) 



and this integral can be explicitly performed in terms of Li 3 , Li 2 and logarithmic functions, 
but the result is quite complicated. For numerical calculations, the integral representation 
of Eq. ATS]) with the adapted quadrature [13| is used to achieve about 64 digits accuracy. 
Having h(n±, 715, n 6 ), the function /(— 2, 0, 0, n 4 , n 5 , uq) is obtained by 

/(-2,0, 0,n 4 , n 5 ,n 6 ) = — [h{n^ n 5 , n 6 ) + n 4 /(-2, 0, 0, n A - l,n 5 ,n 6 )], (79) 



which completes the evaluation of singular Hylleraas integrals. 



VII. SUMMARY 

We have obtained recursion relations for three-electron Hylleraas integrals involving 1 / rf 
and 1/rfj. The initial multi dimensional master integral is expressed in the form of a one 
dimensional integral, which is suitable for high precision numerical evaluation. A set of 
numerical examples for /(— 2, 0, 0, 0, 0, 0) is presented in Table I and for /(0, 0, 0, —2, 0, 0) in 
Table II. Together with regular and extended Hylleraas integrals they allow for the calcula- 
tion of leading QED effects in three electron atoms or ions. Accurate results for the lithium 
atom have already been obtained by Yan and Drake in [3] and several precise measuements 



have been performed in 



L6|, 0| ■ While we intend to verify their result, we aim also to 



calculate QED effects in the Be + ion for the determination of the nuclear charge radii from 
planned isotope shift measurements. Moreover, we aim to calculate higher order QED effects 
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for energy levels along with fine and hyperfine splittings [18( to verify certain discrepancies, 
for example the difference between experimental values and theoretical predictions for the 
2S-3S transition frequency in Li [171. Il9l|. 



TABLE I: Values of the master integral /(— 2, 0, 0, 0, 0, 0) for selected w%,W2, and W3 





W2 


w 3 


/(-2, 0,0, 0,0,0) 


4.0 


1.0 


0.5 


-0.712 305 106 830 898 240 034 381 269 753 396 


4.0 


1.0 


1.0 


-0.699 249 073 328 162 683 321 507 663 896 763 


4.0 


1.0 


1.5 


-0.691 399 920 204 079 424 318 079 206 069 445 


4.0 


1.0 


2.0 


-0.680 940 098 270 060 661 183 207 296 253 486 


4.0 


1.0 


2.5 


-0.667 068 654 115 517 811 166 904 215 261 918 


4.0 


1.0 


3.0 


-0.651 691 691 901 107 499 331 195 660 199 574 


4.0 


1.0 


3.5 


-0.640 014 560 205 730 366 159 157 581 017 121 


4.0 


1.0 


4.0 


-0.627 151 433 983 103 414 178 318 921 819 093 


4.0 


1.0 


4.5 


-0.614 807 615 676 494 534 607 896 770 432 226 


4.0 


1.0 


5.0 


-0.602 911 457 174 426 959 108 113 476 921 161 


4.0 


1.0 


5.5 


-0.591 697 977 384 241 861 061 839 402 802 568 
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7/?o 
w 2 


7/?Q 

w <5 
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Appendix: Singular two electron integrals 

The definition of the distribution P(l/r 3 ) assumes that it is integrated with a smooth 
function of r. This is not always the case with Hylleraas integrals. One can not apply 
directly the definition of Eq. (129]) . because the integrand is logarithmic in r for certain 
negative powers of electron coordinates. Therefore, here we assume the following recursive 
definition 



and set the boundary condition for large 7, by the requirement that only powers of In 7 are 
present, but no constant term. Using this definition, which is consistent with Eq. (1291) . one 
obtains 



d_ 



r(ni, n 2 , n 3 , a, p, 7) = T{n 1} n 2 , n 3 + 1, a, (3, 7) 



(80) 



T(n 1 ,n 2 ,n 3 ,a,l3,j) = 




— e 



ar 1 -(3r 2 -fri2 _ni—l r n 2 ~l r n 3 -l 
'l '2 '12 



(81) 



r(o, 0,0,^/3, 7 ) 



1 



(82) 



(a + /3)(a + 7) (P + 7) 



r(0,0,-l,a,/3, 7 ) 



(a-P)(a + /3) V7 + /3/ 



1 1^1±^ 



(83) 
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wnn On,R^\ (/g + 7) ln(/? + 7 )-(a + 7 ) ln(g + 7) 1 
1 (0,0, -2, a, ^,7) = — — h 



(a 2 - (3 2 ) 



a + (3 



(84) 



r(o, -1,-1,0,0,7) = 1^(1 + ^ + ^(1 + ^ 



(85) 



r(0,-l,-l,a,/3,0) = — 

2 a 



--.n(l + -) In^-L^l-- 



a 



Li,(-- 



(86) 



r(0,-l, -l,a,/3, 7) 



J_ - 2 

2a 



6 2 \a; + 7/ V a + fij V a + 7 



(87) 

T(0, -1, -2, 0, /3, 7) = Li 2 (-1) - 1 ln(l + + l - In 2 - ln(/3 + 7) + y + 1 (88) 



r(o, -1,-2,0^,0) 



l Li2 (-f) + l Inaln ( 1 + l) 



13. 



+(^-i)M« + /J) + ^ + i 



(89) 



r(0,-l,-2,a,/3, 7 ) 



I ( , _ 2) Ll2 f 1 - £±2) - 1 fi + 1) LiJ 1 - 

2V a/ V a + /3J 2\ a) V a + 7 

-i( 1 + 3 ln2 (S) + ^ ln > + «- ln ( Q+ « 



(90) 
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